DMRG Numerical Study of the Kagome Antiferromagnet 
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We numerically study the spin-i antiferromagnetic Heisenberg model on the kagome lattice using 
the density-matrix renormalization group (DMRG) method. We find that the ground state is a 
magnetically disordered spin liquid, characterized by an exponential decay of spin-spin correlation 
function in real space and a magnetic structure factor showing system-size independent peaks at 
commersurate antiferromangetic wavevectors. We obtain a spin triplet excitation gap AE(S = 1) = 
0.055 ± 0.005 by extrapolation based on the large size results, and confirm the presence of gapless 
singlet excitations. The physical nature of such an exotic spin liquid is also discussed. 

PACS numbers: 75.10. Jm,75.50.Ee,75.40.Mg 



Novel magnetic properties and the possible existence of 
exotic spin liquid states [lj in low-dimensional spin-i sys- 
tems have attracted intensive attention in recent years. 
It has been established that spins in the ground state of 
the two dimensional (2D) nearest-neighbor Heisenberg 
antiferromagnet (HAF) model 

H = J2 Si • Sj , 

are still ordered on square Q and triangular 0, 0| lattice 
systems. However, spin liquid states are likely to be 
found in some geometrically more frustrated systems [H, 
like the kagome lattice, which may be seen as a diluted 
triangular lattice (see Fig. [T|) with larger geometrical 
frustration and lower coordination number than the tri- 
angular lattice. Earlier exact diagonalization (ED) stud- 
ies 0, 0, HI suggest that the kagome antiferromagnet has a 
short-range spin correlation and a possible finite spin gap 
~ 0.05 when the finite-size results (up to N = 36 sites) 
are extrapolated to the thermodynamic limit. Within 
the spin gap, a large number of singlet excited states are 
also identified. [H, [a| Recently, algebraic vortex liquid and 
Dirac spin liquid with gapless Dirac fermion excitations 
have been also proposed. 1^, 11 1 Such a Dirac spin liquid 
state has a reasonably good variational energy, ll| but 
the vanishing spin triplet gap is in contrast to the ED 
result. While the discrepancy may be attributed [ll| to 
the uncertainty of the finite size effect in the ED, alterna- 
tively a finite spin gap can be also gained in the Dirac spin 
liquid state via an instabilitvfl2| towards a valence bond 
crystal (VBC) state with a broken translational symme- 
try. Earlier on, Zeng and Marston[l]| also proposed that 
the ground state of the kagome HAF appears to be a 
VBC state with a 36-site unit cell, which is supported 
by the scries expansions. [3] So far the precise nature 
of the HAF on the kagome lattice in the long-wavelength 
and low-energy regime remains unsettled. 

Experimentally the newly synthesized Herbertsmithite 
ZnCu3(OH) 6 Ci2 has brought tremendous excitement to 



this field, in which the spin- ^ copper ions form layered 



Ni a! 

FIG. 1: (color online) Sketch of a three-leg kagome lattice 
with total number of sites N — 3 x Ni x N2 and number 
of unit cells Ni x N 2 = 4 x 3. Here ai = (2, 0) and a 2 = 
(1, v3) are two primitive vectors of the unit cell including 
three inequivalent sites (e.g., 1, 2, 3). 



kagome lattices. The absence of the magnetic order- 
ing has been established based on the neutron scatter- 
ing measurement down to 50mK, as compared to a 
relatively high Curie- Weiss temperature (~ 300 A'). The 
magnetic measurements 15, 3, 17, 3] also suggest that 
there is no signature of a finite spin gap seen in the 
experiment, which seems consistent with an algebraic 
spin liquid, but contrary to a short-range spin liquid 
state with a finite triplet gap. However, possible im- 
purity spins outside the kagome layers, caused by substi- 
tutions of nonmagnetic Zn sites with Cu, or the presence 
of Dzyaloshinsky-Moriya (DM) interactions, [ljj may all 
play an important role in order to fully understand the 
experimental results. While the experimental situation is 
still unclear, on the fundamental side, it is highly desir- 
able to reexamine the issues regarding the nature of the 
ground state and low- lying excitations in a pure spin-i 
HAF model on the kagome lattice. 

Due to the geometrical frustration of the kagome lat- 
tice, the quantum Monte Carlo (QMC) method encoun- 
ters the sign problem, whereas the ED calculation is re- 
stricted to small system size. In this Letter, we present 
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a systematic numerical study by employing the DMRG 
method [2p|, with keeping a large number of basis states 
in the DMRG blocks. We find that the ground state is 
indeed a magnetically disordered state, which is charac- 
terized by an exponential decay of the equal-time spin- 
spin correlation function in real space. The correspond- 
ing magnetic structure factor shows small peaks at com- 
mensurate momenta, with near constant peak values in- 
sensitive to the size of the system, in sharp contrast 
to the structure factor of the magnetic ordered state 
on a triangular lattice. Furthermore, we calculate the 
spin triplet gap, which is extrapolated to a finite value 
AE(S = 1) = 0.055 ± 0.05 in the large sample size limit. 
In this spin liquid state, there also exist low-lying sin- 
glet excitations, with their gap approaching zero at large 
sample size limit. Our calculations strongly hint that 
the ground state may be described by a resonating va- 
lence bond (RVB) spin liquid with short-range antiferro- 
magnetic correlations and weak bond-bond and chirality- 
chirality correlations without explicitly breaking transla- 
tional and rotational symmetries. 

We consider a kagomc lattice with finite length vectors 
iViai and as shown in Fig. [TJ Here ai = (2, 0) and 

&2 = (1, y/o) are two primitive vectors of the unit cell 
which includes three lattice sites on a triangle. The total 
number of sites is N = 3 x Ni x N 2 , with the number 
of unit cells Ni x N 2 . We will extend the calculation 
from N = 36 (the maximum size for ED) to much larger 
sizes with different geometries, up to JV = 3 x 16 x 4 
(192 sites), using the DMRG method. To test the per- 
formance of the DMRG method in the 2D spin systems, 
we have compared our results with the ED up to N = 36 
sites for various lattices (including triangular, square, 
and kagomc lattices) and obtained accurate ground state 
energies with errors smaller than 0.01%. For present 
study, we keep up to m = 4096 states in the DMRG 
block for most systems with more than 24 sweeps to get 
a converged result, and the truncation error is of the 
order or less than 10 -5 . We make use of the periodic 
boundary condition (PBC) to reduce the finite-size effect 
for a more reliable extrapolation to the thermodynamic 
limit. 

We first present the DMRG result for a system with 
N = 48 sites {N x = N 2 = 4). In Fig. [^a), we show the 
ground state energy Eq as a function of m — the number 
of states kept in each block (the dimension of the Hilbcrt 
space = Am 2 ). The ground state energy is extrapolated 
to —20.958 at large m limit and the estimated error at 
in = 4096 is about 0.16%. Similarly, the lowest energy E\ 
in the total spin 5=1 sector is also shown in Fig. [Ha). 
Define the spin triplet gap AE(S = 1) = E^S = 1)-E , 
As plotted in the inset of Fig. [UJa), such a spin gap starts 
to saturate at m > 2000 and approaches the value 0.145 
at large m. In Fig. [^b), similar results for a larger 
system with N = 108 (3 x 12 x 3) are also shown, where 
we find a slightly larger spin gap at 0.163 for this larger 



-20.5 



-20.6 



E? 

CD 
LU 



■20.8 



-20.9 
-46.8 





-46.9 - 




-47.0 - 






■ 




CD 


-47.1 - 


£Z 




LU 


-47.2 - 




-47.3 - 




-47.4 - 



* E1 

• E0 



LU 

< 



0.15 
10 



1500 3000 4500' 

m 

* * * * 



(a) N=3x4x4 



LU 0.1 

< 



o O OO o o 



t>oo o 



1500 3000 4500 

m 



(b) N=3x12x3 



1000 



2000 



3000 



4000 



m 



FIG. 2: (color online) The ground state energy Eo (solid cir- 
cles) and the excitation energy E\ (solid stars) in the total 
spin S — 1 sector are shown as a function of m (the number 
of states kept in each block) for the system with JV = 3x4x4 
in (a) and N = 3 x 12 x 3 in (b). The energy gap for spin-1 
excitation gap Ai? = E\ — Eo (open circles) as a function of 
m for these two systems are shown in the insets. 



but narrow system. 

A systematic size dependence of the spin gap is shown 
in Fig. [3) In the main panel, the spin gap AE(S = 1) 
(solid squares) vs. 1/N is plotted with N = 3 x 4 x 3, 
3x6x3, 3x4x4, 3x6x4, 3x8x4, 3x6x5 
and up to 3 x 6 x 6 = 108, together with the results 
of the ED[t| (open circles) at smaller sizes (note that 
TV = 36 site system in the ED has a different geometry 
as compared to the N = 3 x 4 x 3 system in the present 
calculation). All these data follow nicely a straight line 
shown in Fig. [3j which allows us to extrapolate the spin 
gap to a finite value AE(S = 1) = 0.055 ± 0.005 in the 
thermodynamic limit. Note that all the data presented 
in the main panel are for the systems close to square- 
like with the aspect ratio a = N1/N2 in the range of 
1 < a < 2. The corresponding ground state energies per 
site eo and the spin gaps AE(S = 1) for the various 
system sizes at a given m = 4096 are listed in Table I. 
Furthermore, in the inset of Fig. [31 the spin gap vs. 1 /N 
for 3-leg (JV = 3 x iVi x 3) and 4-leg (N = 3 x N r x 4) 
systems with TVi = 4 — 12 (thus including larger a's) are 
also present for comparison. In general, the spin gap of 
the 4-leg systems is smaller than that of the 3-leg systems 
due to the finite-size effect, consistent with the behavior 
shown in the main panel for the more square-like systems. 
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FIG. 3: (color online) The spin gap AE(S = 1) for square-like 
systems (see text) at different system sizes obtained from the 
ED0] (open circles) and DMRG (solid squares) with N = 36- 
108. The singlet excitation gap AE(S = 0) = Si (5 = 0)- E 
is also given (solid stars). Inset shows the spin gap for both 
3-leg and 4-leg systems. 



Besides the spin triplet gap AE(S = 1), the lowest 
singlet excitation energy AE(S = 0) = Ex(S = 0) — E Q 
is also shown in the main panel of Fig. [3] (solid stars), 
whose magnitude is much smaller than AE(S = 1) and 
approaches zero with increasingsample size. This is con- 
sistent with the ED results^, Q at smaller systems, in 
which a large number of singlet states below the spin 
gap, growing with the system size, are identified. In con- 
trast to the finite spin triplet gap, such vanishing singlet 
excitation energy indicates that the low-lying singlet ex- 
citations will play a dominant role in the low-temperature 
thermodynamic properties like specific heat. 

TABLE I: The ground state energy per site eo and spin gap 
AE(S = 1) for the square-like kagome lattice, obtained by the 
DMRG with keeping m — 4096[2l[ basis states in one block. 



N 


eo 


AE(S = 1) 


3x4x3 


-0.43898 


0.188 


3x6x3 


-0.43875 


0.164 


3x8x3 


-0.43867 


0.165 


3 x 10 x 3 


-0.43868 


0.163 


3 x 12 x 3 


-0.43865 


0.163 


3x4x4 


-0.43591 


0.140 


3x6x4 


-0.43564 


0.122 


3x8x4 


-0.43556 


0.112 


3 x 10 x 4 


-0.43552 


0.114 


3x6x6 


-0.43111 


0.105 




FIG. 4: (color online) The spin-spin correlations \{SoSr)\ 
along the ai direction with N = 3 x 10 x 3 and N = 3 x 10 x 4. 
The error bar represents the mean square deviation of all the 
equivalent sites. The straight line is a fitting to an expo- 
nential function \{SoSr)\ = j4exp(— r/r). The system size 
dependence of the correlation length r is shown in the lower 
left inset. The spin-spin correlations (SoSy ) are also given in 
the insets. 



To characterize the ground state, the spin-spin corre- 
lation function \(SqS*)\ is presented in Fig. 0] for two 
systems with N = 3 x 10 x 3 and 3 x 10 x 4, respec- 
tively. Here r is the distance between the two sites along 
the ai direction in units of the lattice constant and the 
error bar denotes the mean square deviation for all the 
equivalent pairs of sites. Fig. 0] shows that the results 
are well fitted by the straight lines representing an expo- 
nential fit: |(5oi5^)| = Aexp(— r/r) with r as the spin 
correlation length whose size is insensitive to the number 
of legs and is about 0.8 lattice spacing for both systems 
(In the left lower inset of FigHJb), r as a function of N is 
shown for a few systems up to N = 144). These results 
clearly illustrate that the ground state is magnetically 
disordered with no long-range correlations. Furthermore, 
(SqS*) itself exhibits short-range antiferromagnetic oscil- 
lations commensurate with the lattice constant along the 
ai direction as shown in the top right insets of Fig. 2] 
Completely similar results are found for the transverse 
spin-spin correlation (S^~S~) due to the spin rotational 
symmetry. 

To further describe the short-range magnetic correla- 
tions, we calculate the static structure factor S z (q) = 
j? Ey e ici < r '- r ^(SfSp and present the results in Fig. E 
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FIG. 5: (color online) The static structure factors for the 
kagome systems with 7V = 3x4x4in(a);3x6x6in (b); 
The peak values S z (n,n) vs. 1/N are plotted in (c); In (d), 
the structure factor S z (q) for a triangular HAF N — 6 x 6 
system is also shown for comparison. 



(a) and (b) for system sizes JV = 3x4x4 and TV = 3 x 
6x6, respectively, where q are allowed magnetic wavevec- 
tors with components (91,92) = 27r(ni/TVi, 712/TV2) along 
the directions of two primitive basis vectors in the re- 
ciprocal lattice. These figures show that S z (q) exhibits 
small peaks at q* = (tt, 0), (0, w) and (n, tt), indicating 
the dominant short-range antiferromagnetic correlations 
between the nearest neighbor sites. It is important to 
observe that the peaks remain at the same small value 
~ 0.44 without changing much with increasing system 
size, as clearly illustrated by Fig. [5] (c) for the peak 
values at q* = (tt, 7r). Such weak and size-independent 
peaks are in sharp contrast to the structure factor of a 
magnetic ordered system in the triangular HAF model as 
shown in Fig. [5] (d) with TV = 6 x 6 sites, where sharp 
peaks appear at q* = (4p ^) and 4^). 

Finally, we have also checked the bond-bond and 
chirality-chirality correlations and found that both are 
short-ranged with exponential-decay behavior. Thus, 
there seems no explicitly broken translational or rota- 
tional symmetry to be responsible for the gapless sin- 
glet excitation found in such a system. However, for 
our DMRG calculation with finite system size (up to 
192 sites), it is difficult to detect the possible VBC or- 
dering with extremely large unit cell of 36-sites.[H, [3| 
On the other hand, the overall features of the structure 
factor in Fig . [5] (a) and (b) are quite similar to those 
calculated 22j based on the Gutzwiller projected Dirac 
spin liquid state, [ill] although the spin gap vanishes in 



the latter. One may thus conjecture the ground state for 
the kagome HAF be described by an RVB state, which is 
similar to the projected Dirac spin liquid state at short 
ranges. But at long ranges it will have a finite spin triplet 
gap because of the finite size of spin RVB pairing, and 
the gapless singlet excitations are Goldstone modes orig- 
inated from the broken U(l) gauge symmetry due to the 
RVB condensation, like in a charge-neutral superconduc- 
tor. We shall use the variational QMC method to further 
study such kind of spin liquid states and compare with 
the numerical results elsewhere. 

In summary, we have numerically studied the ground 
state properties and low-lying excitations of the kagome 
antifcrromagnct using the DMRG method. Our results 
provide strong evidence that the ground state is a spin 
liquid with only short-range antiferromagnetic correla- 
tions without a magnetic order or other translational 
or rotational symmetry breaking. The spin triplet ex- 
citation has a gap extrapolated to a finite value in the 
thermodynamic limit, but the singlet excitation remains 
gapless. The nature of such a spin liquid state has been 
discussed based on the numerical results. 
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